PPLICATION NO. 09/826,118 
INVENTORS: Urbain Alfred von der Embse 

RECEIVED 

MAY 2 8 2004 

Technology Center 2100 



TITLE OF THE INVENTION 



w Wavelet Multi-Resolution Waveforms 



APPLICATION NO. 09/826,118 

INVENTION: Wavelet Multi-Resolution Waveforms 

INVENTORS: Urbain Alfred von der Embse 



STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 



Not Applicable. 



APPLICATION NO. 09/826,118 

INVENTION: J^frft Wavelet Multi-Resolution Waveforms 

INVENTORS: Urbain Alfred von der Embse 



SEQUENCE LISTING ^ RECEIVED 

MAY 2 8 2004 

Not Applicable. 

Technology Center 2100 



APPLICATION NO. 09/826,118 

INVENTION: Wavelet Multi-Resolution Waveforms 

INVENTORS: Urbain Alfred von der Embse 



INCORPORATION-BY-REFERENCE OF MATERIAL 
SUBMITTED ON A COMPACT DISC 



Not Applicable. 



APPLICATION NO. 09/826,118 

TITLE OF INVENTION: ^iCw Wavelet Multi-Resolution Waveforms 
INVENTOR: Urbain Alfred von der Embse 



5 

BACKGROUND OF THE INVENTION 

I. Field of the Invention 

10 

The present invention relates to CDMA (Code Division 
Multiple Access) cellular telephone and wireless data 
communications with data rates up to multiple Tl (1.544 Mbps) and 
higher (>100 Mbps), and to optical CDMA. Applications are 
15 mobile, point-to-point and satellite communication networks, data 
compression, pattern recognition, media image compression and 
processing, and radar More particularly, the present invention 
relates to Wavelet waveforms and filters for multi-resolution 
applications . 

20 

II. Description of the Related Art 

Wavelets are waveforms of finite extent in time (t) and 
25 frequency (f) over the time-frequency (t-f) space, with multi- 
resolution, scaling, and translation properties. Wavelets over 
the analog and digital t-f spaces respectively are defined by 
equations (1) and (2) . 

30 Continuous wavelet 

, . i i-i/2 J — b . - . 

a 



1 



25 



Digital or discrete wavelet 

Vab( n ) = \ a \ H ) ( 2 ) 

a 



where the two index parameters "a,b" are the Wavelet dilation and 
5 translation respectively or equivalently are the scale and shift, 
"t" is the continuous time variable, and "n" is the discrete 
sample variable for the digital representation assuming that the 
Wavelet is sampled at "T" second intervals. The i|/ "mother" 
wavelet \\f is a real and symmetric localized function in the t-f 
10 space at baseband (dc, f=0) . and generates the doubly indexed 

i r 1/2 

Wavelet v|> ab where the scale factor keeps the Wavelet norm 

invariant under parameter changes xx a,b", and norm is the square- 
root of the energy of the Wavelet response. 

15 Wavelets in digital t-f space are an orthogonal basis with 

the choice of the parameters "a,b" equal to a=2" p , b=qM2 p where 
y p,q" are the new scale and translation parameters and W M" is the 
spacing or repetition interval T S =MT of the Wavelets at the same 
scale "p". Wavelets at "p,q" are related to the mother Wavelet 

20 by the equation 

¥ M (n) =2-^ V (2- p n-qM) (3) 

Wavelets satisfy the orthogonality equation with a correlation 
value equal to "1" when suitably scaled: 



IV P , q Vk,m =1 iff both p = kandq = m 

n 

= 0 otherwise (4) 



FIG. 1 is an N point t-f space which is an N sample 
window of a uniform stream of digital samples at the rate 1/T 
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Hz (1 /second) corresponding to a "T" second sample interval 
and illustrates a Wavelet representation or "tiling" with 
Wavelets that are designed analytically or by an iterated filter 
construction. The t-f space is covered or tiled by a set of 
5 Wavelet subspaces typiP = 0,l v ..,/w-l} 1 where N=2 m . Each Wavelet 

subspace W p 1 at scale "p" consists of the set of Wavelet time 
translations {q = 0,1,..., N/2 p+i -l} over this subspace. These 
Wavelet subspaces are mutually orthogonal and the Wavelets within 
each subspace are mutually orthogonal with respect to the time 
10 translates. This N-point t-f space extends over the time 
interval from 0 to (N-l)T 2 where T is the digital sampling 
interval, and over the frequency interval from 0 to (N-l) 3 in 
units of the normalized frequency fNT 4. 

15 FIG. 2 is an iterated filter bank used to generate the 

Wavelets which cover the t-f space in FIG. 1. Each filter stage 
5 consists of a high pass filter (HPF) and a low pass filter 
(LPF) . Output 6 of the LPF is subsampled by 2 which is 
equivalent to decimation by 2. This t-f space space is an N- 

20 dimensional complex vector metric space V 7. At stage m in the 
iterated filter bank, the remaining t-f space V m 8 is 

partitioned into V m+1 9 and the Wavelet subspace W m+ i 10. 
Scaling functions and Wavelets at each stage of this filter bank 
satisfy the equations 

25 

LPF p : ^(n)=2 1/2 Ih q ^(2n-q) Vq (5) 

q 

HPF p : v|/(n)=2- ,/2 Zg q ^(n-q) Vq (6) 

q 

where 9 is the scaling function, \\f is the Wavelet, HPF P 
coefficients are { h 9 Vq } , LPF P coefficients are {g^, Vq}, and 
30 the equations apply to the stages 0,l,...,m-l. The scaling 
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function (p(n-q) in the Wavelet equation vj/(n) can be changed to 
cp(2n-q) to reflect a downs ampling 12 by the factor 2 if it is 
desirable to reduce the sample rate to the minimum required for 
the subspace and whereupon the 12 would also appear in FIG. 2 as 
5 the final stage of signal processing for each of the Wavelet 
subspaces W 0 ,W lf . . . ,W m 

HPFp and LPF P are quadrature mirror filters (QMF) with 
perfect reconstruction for our application. This means they 
10 cover the subspace V p with flat responses over the subband 
frequency including the edges of the frequency subband, and the 
HPFp coefficients are the frequency translated coefficients for 
the LPF P : {g q =(-l) q h q/ Vq}. 

15 Wavelet design using iterated filter bank starts with the 

selection of the scaling functions. Starting with a primative 
scaling function, one can use the iterated filter construction 
given by equations (5) and (6) to derive successive 
approximations to a desired scaling function which has properties 

20 that have been designed into it by the selection of the filter 

coefficients {h q ,Vq} at each level of iteration. The Wavelets 

can be derived from these scaling functions using the iterated 
filter construction in FIG. 2 or scaling equations (5) and (6) . 
Another use of the iterated filter construction is to design the 
25 scaling functions as Wavelets thereupon ending up with a larger 
set of Wavelets for multi-resolution analysis and synthesis. 

SUMMARY OF THE INVENTION 

30 

This invention is a new and novel design of multi- 
resolution Wavelet waveforms in the frequency domain with the 
capability to be designed to meet specific performance 
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requirements and with a property which allows a single Wavelet 
design to be used for multiple scales. 

Wavelets for linear applications to communications, radar, 
5 data compression, pattern recognition, and media image 
compression and processing, in this invention disclosure are 
assumed to use 1) LS metrics for the passband and stopband 
frequency requirements, 2) quadrature mirror filtering (QMF) LS 
error metric, 3) intersymbol interference (ISI) and adjacent 

10 channel interference (ACI) non-linear LS error metrics, 2) a 
subset of the frequency harmonics as the design coordinates and 
to convert the time domain specification of the LS metrics into 
the frequency harmonics when necessary, and to 3) design the 
Wavelet multi-resolution waveform for no excess bandwidth a=0 

15 where a is the excess bandwidth parameter. These metrics are 
weighed and summed to give a cost function J and iterative 
solution techniques are used to find the set of design harmonics 
which minimize J. These design harmonics generate the mother 
Wavelet and the mother Wavelet is scaled to generate the Wavelets 

20 at different scales specified by the dilation, time translation, 
and frequency translation parameters. 

Non-linear applications of the Wavelets in this invention 
disclosure could require additional metrics and constraints not 

25 included in the linear waveform and filter applications examples 
and are within the scope of this invention disclosure . Examples 
given in this invention disclosure are the constant amplitude 
minimum-shift-keying MSK bandwidth efficient BEM Wavelet 
waveforms for communications and the zero-sidelobe synthetic 

30 aperature radar SAR Wavelet waveforms for radar. 

Scope of this invention disclosure covers decisioning rules 
and metrics other than LS techniques, media (spatial dimension) - 
f space, pattern recognition features-f space, satellite and 
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cellular t-f-beam space, and optical and laser transmission 
time-wavelength space . 

5 BRIEF DESCRIPTION OF THE DRAWINGS AND 

THE PERFORMANCE DATA 

The above-mentioned and other features, objects, design 
algorithms, and performance advantages of the present invention 
10 will become more apparent from the detailed description set forth 
below when taken in conjunction with the drawings and performance 
data wherein like reference characters and numerals denote like 
elements, and in which: 

15 FIG. 1 is an N-point t-f space extending over the time 

interval (0, (N-l)T] and the frequency interval (0, (N-l) /NT] , 
which is tiled or covered by a set of orthonormal Wavelets at the 
scales p=0, l,...,m-l, 

20 FIG. 2 is a Wavelet iterated filter bank used to generate 

the set of Wavelets which tile or cover the t-f space in FIG. 1, 

FIG. 3 is the power spectral density (PSD) template for the 
the stopband and passband used to construct the LS error metrics, 

25 

FIG. 4 is the flow diagram of the LS error metrics leading 
to the construction of the LS error cost function J equal to the 
sum of the weighted LS error metrics and which is minimized 
through selection of the design harmonics to find the optimal 
30 mother Wavelet, 

FIG. 5 is the Matlab 5.0 code for the LS recursive solution 
algorithm used to design the Wavelet waveform in FIG. 6 using an 
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eigenvalue technique, lists the Wavelet time response and the 
design harmonic values, and gives the code for scaling the 
design to a multi-resolution application specified by dilation 
(scale), time translation, and frequency translation 

5 parameters, 

FIG. 6 plots the dc PSD in dB units for the frequency 
response of the mother Wavelet at dc and the square-root raised- 
cosine (sq-rt r-c) waveforms with excess bandwidth parameter a= 
10 0.22, 0.40 vs. the normalized frequency in units of symbol rate 
1/T S . 

FIG. 7 plots the dc PSD in dB units for the frequency 
response of the mother Wavelet at dc and an optimized Gaussian 
15 minimum shift keying (GMSK) waveform vs. the normalized 
frequency fT s . The Wavelet is designed using the LS algorithms 
with additional constraints that reflect the application to the 
minimum shift keying constant amplitude bandwidth efficient 
modulation (BEM) waveform, 

20 

FIG. 8 plots the amplitude of the dc radar ambiguity 
function for a synthetic array radar (SAR) application for the 
mother Wavelet at dc and an optimized unweighted chirp waveform 
vs. the normalized frequency fT p and normalized time t/T c where 
25 the pulse time T p and chip time T c are both equal to the symbol 
interval T s . The Wavelet is designed using the LS algorithms 
with the constraint that . the SAR waveform have no sidelobes. 

30 



35 
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DISCLOSURE OF THE INVENTION 



10 



15 



Wavelet multi-resolution waveforms incorporate a frequency 
translation, are complex, are designed in frequency harmonic 
coordinates, a single mother Wavelet design is used for all 
scales, and in the t-f space are realized as finite impulse 
response (FIR) filters and waveforms. 



1. Wavelet multi-resolution waveform 

Multi-resolution Wavelets at the scales "p,q,r" are related 
to the mother Wavelet at dc by the equation 



\|/ pqr (n) =2-^i|/(2- p n-qM) exp [- j2nf c (p, r ) (2" p n) (2 P T) ] (7) 



where f c (p,r) is the center frequency of the frequency translated 
mother Wavelet at dc as a function of the scale "p" and 

20 frequency index M r" and the phase factor (2 P T) reflects the 
observation that the sampling interval W T" is increased to "2 P T" 
under the scale change from "p=0" to "p=p" The purpose of the 
frequency index "r" is to identify the center frequencies of the 
Wavelets at the scale "p" and translation "q" in the t-f space. 

25 These Wavelets satisfy the complex orthogonality equations 



JW V Uv =1 iff P = k,q = m,r = v (8) 



= 0 otherwise 



where is conjugation. Wavelet basis vectors for the metric 

30 F=t-f space consist of a subset of the admissible set of scaled 
and translated Wavelets { y M , r , Vp,q,r } derived from the mother 

8 



Wavelet \|/ as per equation (7) . An admissible Wavelet waveform 
is any combination that covers F=t-f space. The discrete 
Fourier transform (DFT) harmonic coefficients- {\|/ k } and the 
continuous Fourier transform spectrum i|/(<o) are defined by the 
5 equations 

y/(co) = I>(n) e" ia)n Fourier Transform (9) 

n 

¥k - £?Kn) DFT harmonic coefficients 

n 

where 

y,(a>) = (1/N')2> k Se^'* 

k n 

= 2> k sin((^/2-^k/N*)N*)/ N , sin(^/2-^k/N f ) 

k 

= X^k [Harmonic interpolation for "k"] 

k 

Tirto i2«kn/N t 

rr N , — e 

{k} - DFT frequency or harmonic coefficients such that f k — k/N'T 

where f k is the harmonic frequency corresponding to "k" 
N f = length of i//(n) 



10 Orthonality and no excess bandwidth <x=0 are properties 

which are asymptoticly approached by our Wavelets to within 
design accuracies inherent in communications and radar. The 
Wavelets derived from the mother Wavelet are designed to be 
orthogonal over both time translates "MT" and frequency 

15 translates "l/LT" which respectively correspond to the Wavelet 
symbol spacing T S =MT and the adjacent channel spacing 1/T S . This 
means the orthogonal spacing of the waveforms in V are at the 
time-frequency increments (MT, 1/LT) = (T s . 1/T S ) . In the 
interests of constructing our multi-resolution Wavelets to cover 

20 V it will be convenient to assume that M,L are powers of 2. We 
need the following definitions for the parameters and 
coordinates . 
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Parameters and Coordinates (10) 
N* = Length of \|/ which is an even function about the center and 
which spans an odd number of points or samples 
= ML + 1 where M, L are assumed to be even functions for 

convenience of this analysis 
= Number of points of v|/ 
M = Sampling interval f or 

= Spacing of y/ for orthogonality 
L = Length of v|/ in units of the sample interval M 

= Stretching of y over L sample intervals 
n = n 0 + 

= partitioning into an index n 0 over the sample 

length n 0 = 0,1, ...,M- land an index n, =0,1...,L-1 
over the sample intervals 
k =k 0 + k,L 

= partitioning into an index k 0 over the harmonic 
frequencies k 0 =0,1,...,L-1 corresponding to the 
stretching and an index k x = 0,1 ,...,M - 1 over the harmonic 
frequencies corresponding to the admissible frequency 
slots for y/ 

The harmonic design coordinates are the subset of harmonics 
which span the spectral containment. For the example algorithm 

10 development in this invention disclosure it is assumed the 
waveforms are spectrally contained in the frequency interval 1/MT 
corresponding to the filter spacing which means the L harmonics 
{k 0 =0, 1, . . . ,L-1} covering this spacing are the design harmonics. 
For other applications as will be demonstrated later, the 

15 spectral containment is spread out and one must increase the 

10 



subset of design harmonics. The DFT equations for the mother 
Wavelet in (9) when rewritten in terms of the L harmonic design 

coordinates { y/ k , Vk 0 } become: 



5 DFT equations for the mother Wavelet at dc (11) 

V*o = X^ n ) design harmonics 

n 

^(n) =(l/N')I> k W^ n Wavelet 

10 2. Qrthogonal basis for t-f space 

The Wavelet orthogonality property and tiling of the t-f 
space with a=0 will be proven using the theorems of Karhunen- 
Loeve and Mercer. We start by considering the expansion of a 

15 random complex sequence {z(n), Vn} in a series of Wavelet 
coordinates consisting of time translates of the mother Wavelet 
y/ . The sequence {z(n), Vn} is a zero-mean stationary random 
process which is orthonormal over the sample interval "M" and has 
a frequency spectrum which is flat and extends over the frequency 

20 range 1/MT which is centered at baseband corresponding to a zero 
frequency. This means the {z(n), Vn} cover the subspace of V 
corresponding to the scale of our mother Wavelet i|/ and its time 
translates { y/(n -qM) = y/ q (n\ Vq } . In addition, the V is now 

considered to be extended over a time interval which is 
25 relatively large compared to the N-dimensional t-f space in FIG. 
1,2 to avoid end-effects on the analysis. We start by 
approximating the sequence {z(n), Vn} by the {z(«),Vn}) as per 
the following equations. 
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z(n)= £Z„Kn-qM) 

q 

= EZ q V(n 0 + (n,-q)M) 

q 

where the complex coefficien t s \Z n ] are 

l qJ (12) 

derived from the original sequence 
Z q = Zz(n) ^*(n-qM) 

n 

= Ez(n)<f q (n)* 

D 

wherein the the coefficients {Z q/ Vq} are orthonormal as proven 
by the equations: 

Z q Z q .* = Iz(An-qM)z*(An-q'M)^ q (AnV q . *( An ) (13) 

An 

= ^ qq - E |^ q ( n )| 2 since the sequence {z} is orthonormal 

n 

over the sampling interval "M" 
= L with normalization of the energy of y/ 



Equations (12) and (13) together prove the Karhunen-Loeve 
theorem which proves the following equation for the accuracy in 
approximating the stochastic sequence {z(n), Vn} by {z(«),Vn}. 
This accuracy is expressed by the expected xx E(o)" squared error 
"(o)" in this approximation: 

E{z(n)-z(n)\ 2 }= 1 - J>(«o +n,M)\ 2 (14) 

n | 

We need to prove that the right hand side of this equation is 
zero which then proves that the approximating sequence is equal 
to the original sequence in the mean-square sense. In turn this 
proves that the new waveform coordinates { ys 9 Vq } are an 

orthogonal basis for the original sequence {z(n), Vn} with a=0 
which is our goal. 
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The right hand side of equation (14) when set equal to zero 
expresses Mercer's theorem so our goal is to prove Mercer's 
theorem. To do this we use the DFT of y/ in equation (11) and the 
coordinates in (10) to evaluate the right hand side of equation 
(14) . We find 

i - £h» h-^m)! 2 = l-XIW w*-*r 

n | Icq Ar 0 r 

where T = (VLj^W*^ 

= sin(/r Ak 0 )/Isin(;r Ak 0 1 L) 
= 1 for Ak 0 =0 (15) 
= 0 otherwise 
This proves that 

1- Y\ ¥ (n + Wl M)| 2 =0 Vn 0 

which proves that equation (18) reduces to 

E[z(n)-z(n)\ 2 \= 0 (16) 

which as per the above proves that the set of Wavelets is an 
orthogonal basis and from the construction that a=0. This proof 
easily generalizes to the multi-resolution Wavelets at all of the 
scales {p} and time translates {q} and frequency translates {r}, 
and to the expansion of the harmonic design coordinates as 
required by the application. 

3. Single Wavelet design for multi-scales 

We will demonstrate that a single Wavelet design can be 
used for all scales by 1) implementing the multi-scale 
transformations of the mother Wavelet with the Fourier domain 
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design harmonics, 2) by proving that the mother Wavelet design 
remains invariant under scale changes. 



First consider the multi-scale transformation which derives 
5 the Wavelets at the scale and shift parameters "p,q" for the 
mother Wavelet at scale "p=0" and centered at the origin "q=0". 
We begin by extending the parameters and coordinates in equation 
(10), to include both scaling and subsampling or decimation in a 
form that is equivalent to the iterated filter bank construction 
10 which is used to derive the Wavelet using the filter scaling 
functions. Starting with the coordinates at scale "p=0" the 
parameters and coordinates at scale "p=p" are given by the 
equations: 

15 P=0 parameters and coordinates (17) 

n = n 0 + ni M 
where 

n 0 = a 0 + ai 2 + ... + a m -i 2 m ~ 1 
= M points 
20 M = 2 m 

ni = b 0 + bi 2 + ... 

= points spaced at M sample intervals 
p=p parameters and coordinates 
2" p n = n scaled by "2~ p " 

25 = n(l2 p ), n decimated by 2 P :1 

= n 0 (p) + ni M 

where 

n 0 (p) = n 0 scaled by "2~ p " 
= n 0 (l2 p ) 

30 = a p + ap-i2 + ... + a m - p -i2 m - p - 1 

= M points spaced at 2 P sample intervals 

Combining these equations with the analytical formulation in (7) 
and the Fourier domain representation in (11) enables the 
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Wavelets for "p,q,r" to be written as a function of the Fourier 
domain harmonic design coordinates: 



ip w =(2- p/2 /N')Iv|/ ko W, 



rk 0 (n(p)-qM) 
N 1 



i2^f c (p,r)n(p)2PT 



(18) 



5 for all admissible scale, translation, and frequency index 
parameters, "p,q,r". 

Next we need to demonstrate that the frequency domain 
design in (11) remains invariant for all parameter changes and in 

10 particular for all scale changes. This multi-scale property 
expresses the accordion behavior of the design in that the 
Wavelets at different scales are simply the stretched and 
compressed versions of the mother Wavelet with the appropriate 
frequency translation indices. This multi-scale invariancy means 

15 that the design for a M=16 channel filter bank remains the same 
for M=100 or M=10,000 channel filter banks, when the overlap L 
and the performance goals remain constant. To demonstrate this 
invariant property across scales, we consider the multi- 
resolution Wavelet at scale "p" with the other parameters set 

20 equal to zero for convenience "q=0, r=0" and without loss of 

generality. The Fourier domain frequency response \|/(f) can be 
evaluated starting with the original formulation in equation (7) : 
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DFT at "p, q=0, r=0 



(19) 





sin(/r (fN'2 p T-k 0 )) 



N'sin(;r (fN'2 p T-k 0 )/N*) 



= Z^ko [Harmonic interpolation for" k 0 "] 



This only differs from the harmonic representation in equation 
(9) in the restriction of the design coordinates to the subset of 
harmonic coefficients { y/ k ,Vk 0 } and the stretching of the time 

5 interval to "2 P T" corresponding to the scale "p=p". The harmonic 
interpolation functions are observed to remain invariant over 
scale changes upon observing that the frequency scales as 
"f~l/2 p T" which means the ( f requency*time) product remains 
invariant with scale changes as per the fundamental property of 
10 the Wavelets. This means the Wavelet design is an invariant 
across the Wavelet scales which means we only need a single 
design for all scales of interest. 

15 4. Introduction to LS design 

Design of the Wavelets for a uniform polyphase filter bank 
is described in detail in this invention disclosure. The t-f 
space is spanned by the uniform polyphase filter bank consisting 

20 of M channels at the frequency spacing f w =l/MT where T is the 
digital sampling interval, and the Wavelet filter finite impulse 
response (FIR) time response is stretched over L sampling time 
intervals T s . This polyphase filter bank is ideally decimated 
which means the filter ouput sample rate 1/T S is equal to the 

25 channel-to-channel spacing T S =MT which is equivalent to stating 
that there is no excess bandwidth <x=0 . Our design for this 
topology is immediately applicable to an arbitrary set of multi- 
resolution Wavelet filters through the scaling equation (18) 
which gives the design of our Wavelet at arbitrary scales in 

30 terms of our design of the mother Wavelet at dc. 

LS design algorithms are the eigenvalue and the gradient 
search implementation of an LS algorithm and which respectively 
can be reduced to algorithms that are equivalent to the original 

16 



eigenvalue and Remez-exchange Wavelet design algorithms for 
application to a uniform filter bank. Design of the mother 
Wavelet for this polyphase filter bank uses 5 LS error metrics 
consisting of the 2 prior art passband and stopand metrics which 
5 are upgraded for our mother Wavelet design, and the 3 new 
metrics consisting of the ISI, ACI, and QMF, and solve the LS 
minimization problem using as design coordinates the subset of 
harmonic coordinates. 

10 Frequency domain design coordinates are related to the 

Wavelet time domain digital samples or coordinates as follows. 

Mapping of frequency to time (20) 

15 Time domain design coordinates {^(«),Vn} are real and 

symmetric and can be represented by the reduced set {h t (n), 
n=0,l,...,ML/2} 

h t (n) = if/(0) for n=0 

- 2y/{n) for n=l, 2, ...,ML/2 
20 = time domain design coordinates 

Frequency domain harmonic design coordinates {^,Vk 0 } are 

real and symmetric and can be represented by the reduced 
set {h f (k) , k=0,l,,...,L-l} 

h f (k) = y/ kQ for k=k 0 =0 

25 = 2\f/ K for k=k 0 =l,2,...,L-l 

= frequency domain design coordinates 

Mapping of the frequency coordinates {h f (k), k=0, 1 , L-l } 
into the time coordinates {h t (n), n=0, l,...,ML/2} is defined 
30 by the matrix transformation 

h t = B h f 

where 

h f = (h f (0), h f (ML/2) ) fc transpose of column vector 
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h t = (h t (0), ...,h t (ML/2) ) t transpose of column vector 
B = (ML/2 +1) x L matrix 

- [ Bkn ] matrix of row k and column n elements B kn 
B kn = 1 /ML for n=l 

= 2 cos(27ikn/ ML) otherwise 

wherein the N'=ML+1 has been replaced by ML since a single end 
point has been added to the the FIR to make it symmetrical for 
ease of implementation for the example Wavelet being considered 
with a sample at the mid-point which makes the number of samples 
N' an odd number. 



5. Passband and stopband LS error metrics 

Passband and stopband LS error metrics and cost functions 
for FIR waveforms and filtering are derived with the aid of FIG. 
3 which defines the power spectral density (PSD) parameters of 
interest for the passband and stopband of the PSD ^(o)) for 
communications applications. Requirements for radar applications 
include these listed for communications. Referring to FIG. 3 the 
passband 11 of the Wavelet PSD is centered at dc (f=0) since we 
are designing the or mother Wavelet, and extends over the 
frequency range (Op extending from -(Op/2 to +©p/2 12 in units of 
the radian frequency variable co=27cfT 13 where T is the 
digital sampling interval defined in FIG. 1. The frequency space 
extends over the range of f=-l/2T to f=+l/2T which is the 
frequency range in FIG 1 translated by -1/2T so that the mother 
Wavelet is at the center of the frequency band. Quality of the 
PSD over the passband is expressed by the passband ripple 14. 
Stopband 15 starts at the edge 16 of the passbands of the 
adjacent channels +/-co a /2 16 and extends to the edge of the 
frequency band ct)=+/-7r 17 respectively. Stopband attenuation 
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18 at +/-co a /2 measures the PSD isolation between the edge of the 
passband for the mother Wavelet and the start of the passband for 
the adjacent channels centered at +/-(o s 19. Rolloff 20 of the 
stopband is required to mitigate the spillover of the channels 
5 other than the adjacent channels, onto the dc channel. Deadband 
or transition band 21 is the interval between the passbands of 
contiguous channels, and is illustrated in FIG. 3 by the interval 
from (Op/2 to co a /2 beween the dc channel and adjacent channel at 
co a . Wavelet sample rate <g s 22 is the Wavelet repetition rate. 
10 For the LS example algorithms, the Wavelet sample rate is equal 
to the channel-to-channel spacing for zero excess bandwidth. 
Therefore, 1/T S = co s /27cT=l/MT which can be solved to give gv=2tc/M 
for the radian frequency sampling rate of the filter bank which 
is identical to the Wavelet repetition rate. 

15 

We start by rewriting the DFT equations for the dc Wavelet 
in (11) as a function of the {h f (k), k=0, 1,...,L-1} 

if/ {(d) - ^y(n) cos(n<^) 

n 

= c* B h f using (21) and the definition of the vector V (21) 
c = ( 1 , cos(<y), . . cos( (ML/2) co) )' transpose of column vector 

20 

which is equivalent to the equation for y/(o)) in (9) expressed as 
a linear function of the {h f (k), k=0 , 1 , L-l } . However, this 
functional form is more convenient to analyze. A ideal "c" 
25 vector "c r " will be introduced for the passband and the stopband 
in FIG. 3, in order to identify the error residual 5i//(g)) at the 
frequency "<D" in meeting the ideal passband and stopband 
requirements. The ideal PSD is flat and equal to "1" for the 
passband, and equal to "0" for the stopband. We find 

30 

Error residuals for passband and stopband (22) 
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c r = (l,!,...,].)* 1 passband ideal "c" 
= (O e O,... f O) t stopband ideal "c" 
Sc — c r - c error vector 

Sy/(Q)) = <5c l Bh f 

= Residual error in meeting the ideal spectrum at W 

5 

The LS metric for the passband and stopband can now be 
constructed as follows^ for the eigenvalue and the LS optimization 
(or equivalently, the LS algorithm) design algorithms. 

10 Passband and stopband metrics (23) 

J (band) = — - — (|<V(tf>)| Eigenvalue 
band 2 ' 

= h' f R h f Eigenvalue 

= ||<V I LS 

15 where 

band = [0, ©p ) passband 

= (a> s , 7i] stopband 

R = — — |(B t && t B)d6> 
band J d 

= L x L matrix 
20 8v(/ = ( 8\|/((0i) 8i|/(co u ) ) fc 

= vector of error residuals at the 

frequencies ©i , co u across the band 
||(o)|| = norm or length of the vector (o) and which 

includes a cost function for the errors of 
25 the individual components 



where it is observed that the eigenvalue approach requires that 
the LS metrics be given as quadratic forms in the design 



20 



coordinates {h f (k), k=0, 1,...,L-1} whereas with the LS approach it 
is sufficient to give the LS metrics as vector norms with 
imbedded cost functions or an equivalent formulation. 



6. QMF LS error metric 

Wavelet quadrature mirror filtering (QMF) is a scaling 
property that supports iterative filter construction of Wavelets, 

10 tiling of the t-f space, and perfect reconstruction for polyphase 
filtering. QMF metrics express the requirements on the deadband 
that the PSD's from the contiguous channels in FIG. 3 add to 
unity across the deadband [o)p , G) s ] in order that the filters be 
QMF filters. By suitable modification of the error vector 8c, 

15 the previous construction of the passband and stopband metrics 
can be modified to apply to the deadband. We find 

Deadband metrics (24) 

J (deadband) = (|<^K<^)| Eigenvalue 

deadband 

20 (24) continued 

= h'j-R Eigenvalue 

- 1* i 

where 

5c = c r - c(Q)) - c(7i/M - CO) 
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LS 



where c(©)=c as defined in (22) and (23), and c(n/M - co) = c at 
the offset frequency "n/M - (D" corresponding to the overlap of 
the contiguous filters over the deadband. 



30 
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7. ISI and ACI LS orthogonal error metrics 



Wavelet orthogonality metrics measure how close we are able 
5 to designing the set of Wavelets to be orthogonal over the t-f 
space, with the closeness given by the ISI and the ACI. ISI is 
the non-orthogonal error between channel output samples separated 
by multiples of the sampling interval 1/MT seconds where T is 
the sample time and M is the interval of contiguous samples. ACI 

10 is the non-orthogonal error between between channel output 
samples within a channel and the samples in adjacent channels at 
the same sample time and at sample times separated by multiples 
of the sample interval. When observed as noise contributions 
within each sample in a given channel, the ISI is the noise 

15 contribution due to the other received Wavelets at the different 
timing offsets corresponding to multiples of the sampling 
interval. Likewise, the ACI is the noise contribution due to 
the other Wavelets in adjacent channels at the same sampling time 
and at multiples of the sampling interval. 

20 

ISI and ACI errors are fundamentally caused by different 
mechanisms and therefore have separate metrics and weights to 
specify their relative importance to the overall sum of the LS 
metrics. ISI is a measure of the non-orthogonal between the 

25 stream of Wavelets within a channel as per the construction in 
FIG. 3. On the other hand, ACI is a measure of the non- 
orthogonal between the Wavelet within a channel and the other 
Wavelets in adjacent channels. This means the stopband 
performance metric has a significant impact on the ACI due to the 

30 sharp rolloff in frequency of the adjacent channel, and the ACI 
metric is then a measure of the residual non-orthogonal due to 
the inability of the stopband rolloff in frequency from 
completely eliminating the ACI errors. 
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We assume the received Wavelet is identical to the filter 
Wavelets and is transmitted at the filter output sample intervals 
equal to 1/MT seconds. The second assumption means we are 
assuming the receiver is synchronized with the received signal. 
Since there is no information lost by sampling asynchronously 
with the received Wavelet, we are free to make this 
synchronization assumption without loss of generality. ISI 
metrics are derived in the following set of equations. 



H 

Hjcn 



(25) 

\\f(ML/2)) t transpose of column vector 

matrix of elements H kn 
for n=l 
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ISI metrics 
Mapping of h f into \|/ 

V = (y (-ML/2) , ... 
- H h f 
= (ML+1) x L 
= 1 

= 0.5 cos (27Ckn/ML+l) otherwise 
Offset matrix A 

A = L x (ML+1) matrix of elements A kn 

A kn =0 for k=l 



= [0 0 ... y(-ML/2) ... \|/( (L-k)M+l) ] 
12 kM LM+1 



(25) continued 
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ISI error vector 5E 

8E = L x 1 column vector 

= A H h f 
ISI metric 

J (ISI) = 8E fc 8E Eigenvalue 

= Non-linear quadratic function of h f 

= \\SE\\ ls 
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ACI metrics are derived using the ISI metric equations with the 
following modifications. 

ACI metrics (26) 
5 Mapping of h f into v|/ is the same as developed for ISI 

Offset matrix A elements are changed as follows to apply to 
Channel 1: 

A kn = 0 for k=l 

= [0 0 ... \|/(-ML/2) 0^ ... l|/((L-k)M+l) W£~ k) ] 
10 12 kM+1 LM+1 

which means the ACI error vector 8E is 

8E = L x 1 column vector 
= A H h f 

ACI metric for the two contiguous channels 
15 J (ISI) = 2 8E fc 8E Eigenvalue 

= Non-linear quadratic function of h f 

= 2 \\SE\\ LS 

where the factor "2" takes into account there are two contiguous 
20 channels or one on either side of the reference channel 0 in FIG. 
3. Because of the fast rolloff of the frequency spectrum the 
addition of more channels into the ACI metric is not considered 
necessary, although the functional form of the ACI metric in (26) 
allows an obvious extension to any number of adjacent channels 
25 which could contribute to the ACI. 



8 . LS error cost metric J 

30 Cost metric or function J for the LS algorithms is the 

weighted sum of the LS metrics derived in (23), (24), (25), (26). 
The LS algorithms minimize J by selecting the optimal set of 
frequency coordinates {h f (k),Vk} for the selected set of 
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parameters used to specify the characteristics of the dc Wavelet, 
frequency design coordinates, LS metrics, and weights. Cost 
function and design optimization goals are defined by the 
equations : 

Cost function metric J (27) 
J = ^w(metrics) J(metrics) 

JC 

= weighted sum of the LS metrics J (metrics) 

where 

metrics = passband, stopband, deadband, ISI, ACI 
{w (metrics) , \fmetrics} = set of weights 

^w(metrics) = 1 normalization 

JC 

Optimization goal 

Goal: minimize J with respect to the selection of 
the {h f (k),Vk} 



9. LS error design 

Two classes of LS optimization algorithms considered in 
this invention disclosure as examples for solving equations (27) 
are an eigenvalue and a gradient search implementation of an LS 
algorithm. The eigenvalue algorithm uses a non-linear quadratic 
formulations of the LS metrics and the LS algorithm uses the norm 
formulations for the LS metrics. LS optimization is a class of 
search algorithms which find the minimum LS error cost function 
J. 

FIG. 4 is a summary of the LS metrics and the construction 
of the cost function J. Design parameters 23 are the input and 
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output design parameters. Input parameters are the number of 
polyphase channels M or equivalently the number of digital 
samples at spacing T over the symbol interval T S =MT, the length 
of the FIR time response for the Wavelet in units of L which are 
5 the number of digital samples per Wavelet repetition interval T s 
so that the total number of digital samples for the symmetric FIR 
time response is equal to ML+1, number of DFT samples per FIR 
length n_fft for implementation of the LS algorithms, passband 
radian frequency ©p, stopband radian frequency © a , Wavelet 
10 repetition rate in radian frequency oo s/ selection of the set of 
design coordinates {h f } to be used in the optimization, and the 
metric weights {w (metrics) } . Output parameters are the set of 
harmonic design coordinates {h f } that minimize J. Band metrics 

24 are the passband, stopband, and deadband metrics defined in 
15 equations (23), (23), (24) respectively. Interference metrics 

25 are the ISI and ACI metrics defined in equations (25) and 
(26) respectively. LS cost function J 26 is the weighted 
linear sum of the metrics defined for the band 24 and the 
interference 25 as defined in equation (27) . 

20 

10. Matlab 5.0 code for eigenvalue design 

FIG. 5A, 5B, . . . ,5N,50 are the Matlab 5.0 code for the LS 
25 recursive eigenvalue solution algorithm used to design the 
Wavelet Wavelet in FIG. 6. Included are a listing of the frequency 
domain design harmonics, time response, and data plotting 
software used in the optimization, and an example of the scaling 
of this mother Wavelet into a Wavelet with new scaling 
30 (dilation), time translation, frequency translation parameters, 
and the functions used in the code. The Matlab code uses an 
iterative approach to find the optimal frequency domain 
eigenvector solution to the quadratic form for the LS error whose 
error matrix is the weighed sum of the error matrices for the 
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stopband, passband, deadband metrics and the error matrices for 
the ISI and ACI metrics which are functions of the eigenvectors. 
For each iteration the ISI and ACI error matrices are constructed 
with the frequency domain eigenvectors from the previous 
5 iteration. The iteration finds the new frequency domain 
eigenvectors which are then used to update the ISI and ACI 
metrics for the next iteration. This iteration is repeated until 
there is convergence. The Matlab code allows the operator to 
update the metric weights as well as the scenario parameters to 
10 find the optimal choice for the performance goals. 

FIG 5A Matlab code starts with the selection of the design 
parameters in 29 Step 1. Step 1.1 lists the scenario parameters 
which are M=16 digital samples per Wavelet \\f sample interval, 

15 L=16 nominal length of Wavelet vj/ in units of M, N=ML+1 Wavelet i|/ 
length which is represented by N' in equations (10) , fs=l is the 
normalized channel spacing equal to the Wavelet \\f sample rate, 
fp=0.8864 is the passband frequency interval, n_f=16 is the 
number of design frequency harmonics, n_fft=1024 is the number of 

20 digital samples used for calculation of the frequency spectrum 
centered at 0 frequency. ebno=E b /N 0 =6. 0 dB is the design value of 
the energy per bit E b to noise power density N 0 ratio, and 
x_imbal_aci=6 . 0 dB is the assumed channel-to-channel power 
imbalance used in the calculation of the ACI errors. Step 1.2 

25 derives the software parameters from the design parameters in 
Step 1.1. Step 1.3 lists the optimization parameters which are 
the number of iterations n_iteration = 10, and the metric weights 
w(metric) from equation (27) and in FIG. 4 equal to w(pass) = 
alpha_l=l . e-2 for passband, w (stop) =alpha_2=0 . 80 for stopband, 

30 w(ISI)=alpha_3-2.e-3 for ISI, w (ACI) =alpha_4=0 . 5 for ACI, and 
w (dead) =alpha_5=0 for deadband. The number of iterations was 
selected to provide a convergent solution for the weighted LS 
error metrics w (metric) J (metric) in FIG. 5K 38 Step 10 figure (2) 
and for their sum J in figure (1). Weight values were selected to 
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optimize the Wavelet v|/ filter performance in FIG. 6 and in 
figure (3) in FIG. 5L 39 Step 11, Wavelet y ripple, ISI, ACI 
signal-to-noise SNR power ratio losses in figure (4), and Wavelet 
y time response in figure (5). 

5 

FIG. 5B 30 Step 2 are the initialization calculations prior 
to the start of the iteration loop. Step 2.1 calculates the one- 
sided Wavelet V|/ length parameter m=128. Step 2.2 calculates the 
matrix transformation bw_matrix which maps the \|/ frequency domain 

10 design eigenvectors into the one-sided \j/ time response. Step 2.3 
refers to the function pmn in FIG. 5P,5Q 41 in Step 13.1 which 
computes the LS error matrix which is the sampled data integral 
S [8c°8c' ] el© over the band in equation (23) for the the passband 
and stopband metrics J(pass) , J(stop) . The function pmn_d 

15 computes the corresponding integral over the deadband in equation 
(24) and is not listed since it was not used for FIG. 6. Step 
2.4 constructs a structural matrix c_matrix used in the ISI and 
ACI calculations. Step 2.5 constructs the sample rate templates 
for the passband upper edge, sample rate marker, and passband 

20 lower edge for frequency performance evaluation. 

FIG. 5C 31 Step 3 uses the function pmn to compute the LS 
error matrices passband and stopband in Step 3.1, 3.2. Deadband 
LS error matrix deadband is set equal to a null matrix in Step 
25 3.3. Step 3.4 computes the matrix p_total equal to the weighted 
sum of these matrices. Step 3.5 maps this LS error matrix 
p_total in the time domain into the LS error matrix pw_t in the 
frequency domain using the matrix transformation bw_matrix. 

30 FIG. 5C,5D 32 Step 4 initializes matrices used in the 

iteration loop and starts the iteration loop for the iterations 
i_iteration=l, 10 . Step 4.1 calculates the eigenvalue and 
eigenvector which minimize the cost function J in equations (27) 
and FIG. 4 whose LS error matrix is pw_t. Step 4.2 uses the 
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bw_matrix to map the eigenvector in frequency into the Wavelet 
time response i|/=hn and stores the eigenvector as the set of 
Wavelet \\j design harmonics i|/ k =hw_eig (k) for k=0,l,...,15 where i|/ k 
has been defined in equations (9) and in equations (11) upon 
5 recognizing that there are n_f=16=L design harmonics k o =k=0,l,. . 
.,15. There is no reason to use the negative harmonics since the 
spectrum is real and symmetric about the 0 frequency. Step 4.3 
computes the passband peak-to-peak ripple as xripple and the 
stopband attenuation as xstop. 

10 

FIG. 5E 33 Step 5 calculates the weighted passband error 
metric w (pass) J (pass) =beta_pass, weighted stopband error metric 
w(stop) J (stop) =beta_stop, and weighted deadband error metric 
w (dead) J (dead) =beta_dead in equations (27) using the b_vector 
15 format for hn and the LS error matrices passband, stopband, 
deadband. 

FIG. 5E,5F 34 Step 6 calculates the ISI and ACI LS error 
matrices, their metrics, and their error contributations to the 

20 SNR loss. Step 6.1 calculates the ISI LS error matrix w__matrix 
and the corresponding ISI metric J ( ISI ) =errM_isi with Matlab code 
which implements equations (25) , and calculates the ISI residual 
error errV_isi that causes an ISI SNR loss. Step 6.2 calculates 
the ACI LS error matrix w_f_matrix and the corresponding ACI 

25 metric J(ACI) = errM_aci with Matlab code which implements 
equations (26) , and calculates the ACI residual error errV_aci 
that causes a ACI SNR loss. 

FIG. 5G 35 Step 7 calculates the weighted LS error metrics 
30 for ISI and ACI and the LS error matrix pw_t for the next 
iteration. Step 7.1 calculates the weighted ISI LS error metric 
w (ISI ) J (ISI ) =beta_isi and the weighed ACI LS error metric 
w(ACI) J (ACI) =beta_aci in equations (27) and the sum of the 
weighted LS error metrics J=errM_isi. Step 7.2 saves the weighted 
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LS error metrics and their sum J for each iteration. Step 7.3 
updates the pw_t matrix for the next iteration. 

FIG. 5G,5H 36 Step 8 calculates the SNR losses, stores 
5 these losses for each iteration, and completes the iteration 
loop. Step 8.1 computes the SNR loss in dB due to passband 
ripple as xloss_ripple, due to ISI as xloss_isi, due to ACI as 
xloss_aci, and the total SNR loss as xloss_total. Step 8.2 saves 
these SNR losses for each iteration and completes the iteration 
10 loop. 

FIG. 5H,5I,5J,5K 37 Step 9 documents the Wavelet y 
frequency design harmonics and time response generated by the 
Matlab code for the Wavelet in FIG. i|/ 6. Step 9.1 lists the 

15 Wavelet \\? frequency domain harmonic design coordinates 
\|/ k (k)=hw_eig(k) for harmonics k=0, 1,..., 15. Only the positive 
harmonics are listed since the frequency spectrum is real and 
symmetric about k=0 . Step 9.2 lists the Wavelet y time response 
y(n)=hn(n) for n=0,l,2,. . . ,128 with the coordinate frame 

20 centered at the Wavelet peak. Only the positive digital samples 
are listed since the time response is real and symmetric about 
n=0 . 

FIG. 5K 38 Step 10 plots the total weighted LS error metric 
25 J in figured) and the individual weighted LS error metrics in 
figure (2) as functions of the iteration number. After about 6-7 
iterations both the total and the individual weighted LS error 
metrics are observed to stabilize indicating that the selected 
number of iterations n_iteration=10 is sufficient for 

30 convergence. 

FIG. 5L 39 Step 11 plots the Wavelet i|/ frequency response 
in figure (3) which is the same as the frequency response in FIG. 
6, plots the SNR losses due to passband ripple, ISI, ACI and the 
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total loss in figure (4), and plots the time response in figure (5) 
for use in the optimization of the LS metric weights in Step 1.3 
as well as the scenario parameters in Step 1.1. 

5 FIG. 5M 40 Step 12 documents the Matlab code to rescale the 

mother Wavelet \\i in FIG. 6 which was derived in the previous 
Matlab code in Steps 1-11, for new scale (dilation) p, time 
translation q, and frequency translation k parameters compared to 
the values p=0,q=0,k=0 for the mother Wavelet \\t . An example set 
10 of new parameters documented in the code is p=2, q=2, k=3 and 
with the assumption that the digital sampling rate remains fixed 
under the parameter change. A constant digital sample rate keeps 
the frequency band constant under the parameter change. 

15 FIG. 5M 40 Step 12.1 defines the new Wavelet i|/ p , q/k in 

equations (7) ,(18) resulting from these parameter changes, in 
equations (28) 

20 New Wavelet scaled from mother Wavelet (28) 

M> P ,q,k= 2" p/2 ij/(2- p n-qM) exp(j27i k 2" p n /ML) 

Case 1 Scaled (dilated) sampling 2" p n with fixed M 
25 subsamples (decimates) n by the factor 2 P 

• Represent n as 

n = n_0 + 2 P n_p 
n_0 = 0,1, . . ., 2 P -1 
n_p = 0, + /-1, +/-2, . . . 
30 • We observe 2~ p n = n_p 

* This enables the new Wavelet to be written 
H> P , q ,k = 2~ p/2 \|/(n_p-qM) exp(j2rc k n_p /ML) 

Case 2 Sampling n remains fixed and M is 
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rescaleci by the factor 2 P 

• New M is M_new 

M_new = 2 P M 

• This enables the new Wavelet to be written 

5 M/ P , q , k = 2~ p/2 V|/ (n-qM_new) exp(j27i kn/M_newL) 



wherein the new Wavelet equation from (7) ,(18) is re-written 
using the frequency translation parameter k=r and the explicit 
frequency translation expression exp(j2n k 2" p n /ML) which is 
10 recognized as the discrete Fourier transform kernal. 



In Case 1 corresponding to the standard Wavelet assumption, 
the subsampling or decimation of n by the factor 2 A p stretches 
the time between contiguous samples by the factor 2 A p while 

15 keeping constant the number of samples M per Wavelet sample 
interval, and therefore reduces the frequency band by the factor 
2 A p. The new n can be represented as the sum of the two fields 
n_0 = 0,1,. . , (2 A p -1) and n_p = 0, + /-1, +/-2, . . . Scaling 
or equivalently subsampling results in the new n equal to n_new 

20 = 2 A -p n = n_p corresponding to a stretching of the sample 
interval by 2 A p. The new Wavelet length N_new remains constant 
N_new=N using the n_p sampling. 

In Case 2 which is the communications application being 
25 considered, the sampling is held constant corresponding to a 
constant frequency band assumption. The factor 2 A -p is divided 
out in the Wavelet expression which is equivalent to the M being 
scaled to M_new=2 A p M corresponding to increasing the number of 
samples per Wavelet interval by the factor 2 A p. This means Case 
30 2 corresponds to increasing the maximum number of channels by the 
factor 2 A p over the same frequency band. The new Wavelet length 
is equal to N_new=M_newL+l =1025 since M_new=2 A pM=64 . 
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FIG. 5M 40 Step 12.2 derives the new matrix transformation 
bw_matrix_new which maps the \\f frequency domain design 
eigenvectors into the one-sided \y time response. 

5 FIG. 5M, 5N 40 Step 12.3 uses the matrix transformation 

bw_matrix_new to map the design eigenvector hw_eig into the new 
Wavelet in three steps. In the first step, the mother Wavelet is 
translated into the new baseband Wavelet vj/ P/q=0 ,k=o = hn_0 which is 
the mother Wavelet scaled by p and with no time and frequency 
10 translations q=k=0. The second step translate this baseband 
Wavelet by qM_new to generate Vp,q,k=o = hn_l . The last step 
translates this Wavelet v|/ p , q , k=0 = hn_l in frequency to give the new 
Wavelet x|/ p , q , k = hn_new. 

15 FIG. 5N 40 Step 12.4 plots the Wavelet time responses for 

the mother Wavelet v|/=hn and the new Wavelet i|/ p , q , k = hn_new in 
figure (6). Step 12.5 plots the Wavelet frequency responses for 
the mother Wavelet v|/=hn and the new Wavelet i|/ P , q , k = hn_new in 
figure (7) versus the normalized frequency/ (vj/=hn sample rate) and 

20 in figure (8) versus the normalized f requency/ (M/ P , q , k = hn_new 
sample rate) . 

FIG. 50 41 Step 13 documents the functions used in the 
Matlab code. Step 13.1 is the code for the function pmn which 
25 computes the matrix for the J (BAND) in equations (23). FIG. 5P 41 
Step 13.2 lists the code for the function freq_rsp which computes 
the Fourier transform of the input hn versus the normalized 
frequency/ (Wavelet \y=hn sample rate). 
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11. Wavelet applications 

Applications of this new invention to both communications 
and radar include the use of the new Wavelet: 1) for 

5 communications applications which use the square-root raised 
cosine waveform (sq-root rc) which is extensively used for the 
third generation (3G) CDMA communications, 2) for communications 
applications which use the Gaussian minimum shift keying (GMSK) 
waveform for constant amplitude bandwidth efficient (BEM) 
10 applications, and 3) for synthetic aperture radar (SAR) and 
real aperture radar (RAR) applications. 

FIG. 6 illustrates the communications applications to 
filters and waveforms currently using the square-root raised- 

15 cosine (sq-rt rc) waveforms with bandwidth expansion parameter 
a = 0.22 to or =0.40. This notation means that for a =0.22 the 
spectral efficiency is (symbol rate/bandwidth) = 1/1+a =1/1.22 
= 0.82 = 82%. CDMA communications applications for sq-rt rc 
waveforms include current CDMA and 3G CDMA. FIG. 6 plots dc 

20 power spectral density or power spectrum (PSD) in units of dB 42 
versus the frequency offset from dc expressed in units of the 
symbol rate 43 for the Wavelet waveform 44, the sq-rt r-c 
with a=0.22 45, and the sq-rt r-c with a=0.40 46. 

25 

12. Wavelets for minimum shift keying BEM 

Constant amplitude BEM application of the new Wavelet 
30 waveform indicates that it can be used for communications 
applications of GMSK. The GMSK finds applications for 
transmitters which operate their HPA(s) amplifiers in a 
saturation mode in order to maximize their radiated power from 
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the HPA(s), and which require a BEM PSD to avoid excessive 
spreading of the transmitted power. 



FIG, 7 plots the simulation data for the Wavelet waveform 
5 BEM 49 and the GMSK 50 as PSD in dB 47 versus the normalized 
frequency offset from dc expressed in units of the bit rate 48. 
Both Wavelet and GMSK use a a length parameter L=10 where L is 
the length of the phase pulse in units of phase pulse repetition 
rate. 

10 

13. Wavelets for synthetic aperture radar 

Radar RAR and SAR applications of the Wavelet waveform 
15 indicate that it can be used for the chirp waveforms for wideband 
signal transmission, when combined with pseudo-random phase 
codes . 

FIG. 8 plots the results of the simulation for the Wavelet 
20 waveform and an unweighted frequency chirp waveform. Plotted are 
the ambiguity function for the Wavelet waveform 51 and the 
unweighted frequency chirp waveform 52. The dc 2-dimensional 
radar ambiguity function 53 is plotted as a function of the 
frequency offset in units of fT p 54 and the time offset in 
25 units of t/T c 55 where T p is the phase-coded radar pulse length or 
length of the phase code and T c is the phase code chip length. 
The chip length is identical to the waveform repetition interval 
T s so that T c = T s . 

30 
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14. Preferred embodiments 



Preferred embodiments in the previous description is 
provided to enable any person skilled in the art to make or use 
the present invention. The various modifications to these 
embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the 
present invention is not intended to be limited to the 
embodiments shown herein but is to be accorded the wider scope 
consistent with the principles and novel features disclosed 
herein. 
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